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Abstract 

The  liquid  repellence  and  surface  topography  characterstics  of  coatings  comprising  a  sprayed-on 
mixture  of  fluoroalkyl-functional  precipitated  silica  and  a  fluoropolymer  binder  were  examined  using 
contact  and  sliding  angle  analysis,  electron  microscopy,  and  image  analysis  for  determination  of  fractal 
dimensionality.  The  coatings  proved  to  be  an  especially  useful  class  of  liquid  repellent  materials  due  to 
their  combination  of  simple  and  scalable  deposition  process,  low  surface  energy,  and  the  roughness 
characteristics  of  the  aggregates.  These  characteristics  interact  in  a  unique  way  to  prevent  the  build-up 
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of  binder  in  interstitial  regions,  preserving  re-entrant  curvature  across  multiple  length  scales,  thereby 
enabling  a  wide  range  of  liquid  repellency,  including  superoleophobicity.  In  addition,  rather  than 
accumulating  in  the  interstices,  the  binder  becomes  widely  distributed  across  the  surface  of  the 
aggregates,  enabling  a  mechanism  in  which  a  simple  shortage  or  excess  of  binder  controls  the  extent  of 
coating  roughness  at  very  small  length  scales,  thereby  controlling  the  extent  of  liquid  repellence 

1.  Introduction 

A  significant  advance  in  the  field  of  liquid  repellence,  the  creation  of  surfaces  exhibiting 
superoleophobicity  has  generated  much  interest  since  2007. 1  The  unusual  characteristics  of 
superoleophobic  surfaces  in  contact  with  fuels,  oils,  and  greases,  including  contact  angles  in  excess  of 
150°,  droplet  sliding  angles  approaching  zero,  and  robust  metastability  of  the  partially  wetted  (Cassie- 
Baxter)  state,  '  provide  exceptional  liquid  repellence  with  applications  in  machinery,  fabrics,  ’  energy 
efficiency,9  and  protection  from  fouling  and  hazardous  materials.10  After  only  very  sporadic  previous 
reports,  researchers  have  recently  created  many  types  of  superoleophobic  surfaces,  with  new  reports 
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appearing  at  a  very  rapid  pace.  '  ’  ’  ’  ’  '  Fabrication  methods  involve  dip  coating, 
electrospinning,  ’  ’  ’  spray  coating,  and  templating.  Among  these,  spray  coating  provides  the 
advantages  of  great  simplicity,  scalability,  speed,  and  compatibility  with  almost  any  substrate  type  (if 
post-processing  is  limited  to  drying  at  low  temperature  and  appropriate  solvents  are  selected).  Thus, 
there  remains  a  significant  need  for  further  development  of  superoleophobic  surfaces  fonned  by  spray 
coating,  as  well  as  for  simple  means  of  tuning  the  liquid  repellence  of  such  surfaces. 

Herein,  we  describe  the  development  of  a  set  of  materials  and  methods  for  producing  spray- 
coated  surfaces  with  liquid  repellence  characteristics  ranging  from  hydrophobic  to  superoleophobic 
(toward  hydrocarbons  as  light  as  dodecane)  that  address  the  challenges  of  scalable,  affordable 
fabrication  and  straightforward  tuneability.  The  simple  and  robust  process  for  fabricating  these  surfaces 
utilizes  fluoroalkylsilane-treated  precipitated  silica  aggregates  (“FF-silica”),  the  synthesis  of  which  we 
have  previously  reported.  The  coating  formulation  also  includes  a  Viton®  fluoropolymer  binder,  and 
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a  fluorinated  solvent  (AK-225G).  This  combination  resulted  in  a  wide  range  of  stochastic  surface 
topographies  and  very  low  surface  energy  values.  (This  combination  of  materials  has  also  been  widely 
utilized  to  achieve  high  levels  of  abrasion,  wear,  and  scratch  resistance.  )  Roughness  on  multiple 
length  scales  below  10  pm,  along  with  the  propensity  to  form  surfaces  with  re-entrant  curvature,  were 
observed  to  be  key  factors  that  enabled  these  surfaces  to  exhibit  a  wide  range  of  liquid  repellency.  In 
addition,  the  ability  to  control  the  fine-scale  texture  by  altering  the  amount  of  binder  in  the  system  was 
found  to  be  a  convenient  geometric  tuning  mechanism  for  shaping  the  liquid  repellence  characteristics 
of  these  randomly  textured,  sprayed-on  surfaces. 

2.  Materials  and  Methods 

2.1  Preparation  of  Fluoroalkyl-Functional  Silica  (FF-Silica):  Precipitated  silica  aggregates  (HiSil 
233,  from  PPG  Industries)  were  treated  with  (l,l,2,2-tetrahydroperfluorodecyl)-dimethylchlorosilane 
using  a  previously  described  method."  Relevent  physical  properties  of  the  modified  silica  aggregates 
can  be  found  in  Table  S 1  of  Supporting  Information. 

/p\ 

2.2  Spray-Deposition:  Coatings  were  formulated  from  a  stock  mixture  of  5  mg/mL  Viton  ETP-600S 
fluoropolymer  (DuPont)  in  Asahiklin  AK-225G  (l,3-dichloro-l,2,2,3,3-pentafluoropropane;  AGC 
Chemicals  Americas).  FF-silica  was  suspended  into  this  mixture  at  0-90  wt%  silica/total  solids. 
Additional  solvent  was  added  to  maintain  a  total  solids  concentration  of  5  mg/mF.  Samples  were 
prepared  by  spraying  thoroughly  mixed  coating  fonnulations  immediately  onto  1"  silicon  wafers  (P- 
type;  100-orientation,  Wafer  World)  through  an  airbrush  (Paasche,  VFSTPRO)  with  a  1.06  mm 
diameter  tip  using  compressed  air  (25  psi)  .  The  airbrush  was  repeatedly  passed  over  the  substrate 
laterally  at  an  approximate  distance  of  15-20  cm  from  the  substrate  until  20  mF  of  the  coating  mixture 
had  been  deposited.  The  resultant  deposition  level  is  around  is  20  mg/cm".  Following  spray  coating, 
samples  were  air  dried  for  1  hr  followed  by  drying  for  12  hrs  at  60  °C  (ambient  atmosphere)  in  a 
laboratory  oven. 

2.3  Contact  Angle  Measurements:  Advancing  (0adv)  and  receding  ( 0rec)  contact  angles  were  measured 
using  a  VCA  2000  goniometer  (AST,  Inc.)  as  probe  fluid  was  added  or  removed  from  ~  5  pF  sessile 
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droplets.  Sliding  angles  were  measured  with  a  Rame-Hart  Model  590  goniometer  on  -15-20  pL 
droplets.  Measurements  were  made  at  multiple  locations  in  each  sample,  with  the  standard  deviation 
from  all  measurements  being  reported  as  the  characteristic  uncertainty  in  Table  S2  of  Supporting 
Information.  For  comparison,  additional  measurements  on  several  samples  with  different  surface 
composition/probing  liquid  combinations  were  made  on  a  Dataphysics  OCA20  goniometer  equipped 
with  a  TBU90  tilting  stage,  and  no  significant  differences  were  found. 

2.4  Surface  Characterization:  Cross-sectioned  samples  were  prepared  by  embedding  coatings  in 
epoxy  and  then  sectioning  the  embedded  sample  using  a  PT-X  ultramicrotome  (RMC  Products)  using 
Boeckeler,  Inc.  methodology.  Gold-sputtered  cross-section  and  plan  view  (top  down)  samples  were 
imaged  using  an  FEI  Quanta  600  scanning  electron  microscope. 

3.  Results  and  Discussion 

The  wide  range  of  liquid  repellency  of  the  sprayed-on  surfaces  based  on  the  combination  of 
precipitated  silica  aggregates,  fluoropolymer  binder,  and  AK225G  is  illustrated  in  Fig.  1.  At  low  FF- 
silica  loading  levels,  the  Wenzel  (fully  wetted)  state  exists  for  all  liquids  studied,  with  advancing  contact 
angles  that  are  typical  of  highly  fluorinated  surfaces.  ‘  As  the  FF-silica  loading  level  increases, 
however,  the  Cassie-Baxter  state  becomes  prevalent,  starting  at  approx.  20  wt%  FF-silica  for  water  and, 
at  higher  loading  levels,  for  liquids  of  progressively  lower  surface  tension.  As  indicated  by  the  more 
detailed  liquid  repellence  map  in  Fig.  2,  a  generally  simple,  monotonic,  and  gradual  relationship  exists 
between  the  surface  tension  of  the  contacting  liquid  and  the  FF-silica  loading  level  required  for  attaining 
the  Cassie-Baxter  state.  Although  Fig.  2  applies  directly  only  for  the  specific  set  of  spraying  parameters 
reported  herein,  it  should  be  noted  that  the  spray  process  includes  an  annealing  step  above  the  glass 
transition  of  the  polymer  binder,  which  allows  for  relaxation  of  any  residual  stresses  and  some 
equilibration,  even  after  the  solvent  has  evaporated.  As  a  result,  we  expect  that  the  behavior  illustrated 
in  Fig.  2  would  exhibit  a  low  sensitivity  to  variations  in  processing  parameters,  provided  that  a  similar 
overall  density  of  solids  is  deposited  on  the  surface.  The  existence  of  this  simple  and  robust  relationship 
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enables  quantification  of  the  trade-offs  (often  driven  by  the  need  to  simultaneously  satisfy  multiple 


perfonnance  criteria)  required  to  obtain  a  desired  degree  of  liquid  repellence,  and  suggests  that  even 


complex,  randomly  textured  surfaces  may  be  systematically  altered  in  a  straightforward  manner  by 


altering  the  ratio  of  aggregates  to  binder  in  order  to  control  liquid  repellence. 
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Figure  1.  Apparent  advancing,  receding,  and  sliding  contact  angle  measurements  of  water  (y/v  =  72.1 


mN/m),  diiodomethane  (y /v  =  50.8  mN/m),  rapeseed  oil  (y/v  =  35.5  mN/m),  and  hexadecane  (y /v  =  27.5 


mN/m)  on  flat  silicon  wafers  sprayed  with  mixtures  of  fluoropolymer  and  FF-silica. 
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Figure  2.  Map  depicting  the  state  of  liquid  contact  with  FF-silica  surfaces  as  a  function  of  FF-silica 
loading  and  surface  tension  of  the  contacting  liquid;  fdled  circles  represent  Wenzel  states,  shaded 
triangles  represent  Cassie  states  with  contact  angle  hysteresis  >  20°,  unfilled  squares  represent  Cassie 
states  with  contact  angle  hysteresis  <  20°.  The  background  color  is  intended  only  as  an  aid  for 
visualization.  Inset  shows  droplets  of  water  (blue)  and  dodecane  (purple)  on  spray-coated  surfaces 
comprised  of  50  wt%  FF-  silica  (left  image)  and  80  wt%  FF-silica  (light  image)  demonstrating  the 
different  wetting  states  as  a  function  of  silica  loading. 


In  order  to  better  understand  the  basis  for  the  simple  relationship  seen  in  Fig.  2,  the  geometry  of 
the  sprayed-ou  surfaces  was  examined  in  more  detail.  Fig.  3  provides  both  surface  view  and  cross- 
sectional  SEM  images  of  the  surfaces  as  a  function  of  FF-silica  loading,  with  detailed  views  provided  in 
Fig.  4.  An  analysis  of  the  fractal  dimensionality  of  the  surfaces  as  determined  by  the  method  of 
Shibuichi  “  is  shown  in  Fig.  5.  The  images  in  Figs.  3  and  4  demonstrate  that  the  FF-silica  aggregates 
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(which  feature  a  fluoroalkylated  surface  to  ensure  liquid  repellence  and  maintain  good  dispersion  in  the 
spraying  solvent)  retain  sufficient  integrity  to  produce  large  protrusions  from  the  fluoropolymer  layer 
after  deposition  and  film  formation.  It  also  appears  that  the  low  surface  tension  of  the  carrier  solvent, 
perhaps  in  combination  with  the  porosity  afforded  by  the  aggregates,  allows  the  fluoropolymer  coating 
to  become  confonnal,  particularly  as  the  FF-silica  loading  level  is  increased  beyond  approx.  40%.  In 
Fig.  4,  regions  of  conformality  are  indicated  by  filled  arrows,  while  regions  that  lack  conformality  are 
indicated  by  unfilled  arrows.  Since  the  conformality  extends  to  the  base  of  the  protrusions,  it  appears 
the  typical  phenomenon  of  binder  accumulation  in  the  interstitial  junctions  does  not  take  place.  Rather, 
the  protrusions  exhibit  re-entrant  curvature,  a  key  geometric  characteristic  that  enables  the  Cassie  state 
to  be  observed  with  low  surface-tension  liquids.1  In  addition,  the  conformality  transfers  the  multi-scale 
roughness  inherent  in  the  aggregates  to  the  coating  as  a  whole.  Therefore,  the  use  of  silica  aggregates  in 
a  fluoropolymer  binder  with  a  fluorocarbon  carrier  solvent  allows  the  key  characteristic  of  re-entrant 
curvature  to  appear  at  multiple  length  scales  in  sprayed-on  FF-silica  coatings.  In  combination  with  the 
very  low  surface  energy  imparted  by  the  FF-silica  and  the  fluoropolymer,  the  re-entrant  curvature  leads 
to  outstanding  liquid  repellency,  including  superoleophobicity. 
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Figure  3.  Plan  view  (top  panels)  and  corresponding  cross-sections  (lower  panels)  for  FF-silica  / 
fluoroelastomer  composites  at  FF-silica  loadings  of  a)  20  wt%  b)  40  wt%  c)  60  wt%  d)  80  wt%.  Note 
that  the  images  are  all  at  identical  magnification,  so  that  the  scale  bar  applies  to  all  images. 
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Figure  4.  SEM  micro-graphs  prepared  from  FF-silica  /  fluoropolymer  coatings  mounted  in  epoxy  and 
cross-sectioned  (all  at  magnification  3000  x).  Note  the  scoring  of  the  epoxy  due  to  the  sectioning 
process,  which  is  helpful  in  distinguishing  it  from  the  fluoropolymer  binder  that  surrounds  the  silica 
aggregates.  Filled  arrows  indicate  examples  of  regions  where  the  coating  is,  if  present  at  all,  highly 
conformal  to  fine  features  of  the  silica  surface,  whereas  unfilled  arrows  indicate  examples  where  the 
coating  does  not  conform  to  the  fine  features  of  the  silica  surface. 
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Figure  5.  Fractal  dimensionality  of  FF-silica  /  fluoroelastomer  composites  as  a  function  of  FF-silica 
loading 


Although  the  exact  level  of  coating  conformality  is  difficult  to  determine,  it  is  clear  from  the 
images  in  Fig.  4  that  even  features  of  <  1  pm  are  preserved.  In  accordance  with  the  concept  of  a 
minimum  wetted  pore  diameter,  which  strictly  applies  only  to  liquids  for  which  the  equilibrium  contact 
angle  exceeds  90°,  and  which  should  be  on  the  order  of  0.1  pm  for  the  highly  fluorinated  aggregates, 33 
the  surface  of  any  highly  conformal  coating  on  the  aggr  egate  is  also  likely  to  be  only  partially  wet.  In 
terms  of  the  systematic  design  parameter  approach  for  regularly  patterned  surfaces,  a  dimensionless 
parameter  D*  can  be  defined,  based  on  <f>SJOtai,  the  fraction  of  the  liquid-solid-air  composite  interface  that 
is  occluded  by  solid,  as 

D*=  \  I  <j>s, total  (1) 

In  most  calculations  of  D*  involving  discrete  particles  on  a  surface,  it  is  assumed  that  the  surface  of 
each  particle  is  frilly  wetted,  while  the  spaces  between  particles  remain  unwet.  However,  if  the  surface 
is  comprised  of  aggregates,  and,  due  to  fine-scale  roughness,  some  fraction  “w”  of  the  surface  of  each 
aggregate  remains  umvet,  then  D*  can  be  expressed  in  a  modified  form  given  by 
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D*=l/[(l-u)fa0]  =  [1/(1  -u)]Do* 


(2) 


In  Eq.  (2),  the  quantities  with  subscript  0  refer  to  equivalent  values  for  a  surface  comprised  of 
monolithic  particles  with  an  arrangement  identical  to  that  of  the  aggregates.  Because  0  <  u  <  1,  D*  will 
be  as  high  or  higher  for  a  surface  conformal  to  aggregates  as  compared  to  an  identical  surface  that 
conforms  to  monolithic  particles.  The  value  of  u  will  depend  on  specific  details  of  the  surface 
topography  at  the  line  of  contact,  ’4  not  simply  the  porosity  of  the  aggregate,  hence  reliable  a  priori 
methods  for  calculating  u  are  not  readily  available.  For  unifonnly  distributed,  monolithic,  and 
monodisperse  spherical  particles,  however,  the  parameter  Do *  can  be  thought  of  as  the  square  of  an 
effective  spacing  ratio,  that  is  the  distance  between  particles  divided  by  the  particle  diameter  (see 
Supporting  Information,  Section  S3).  Although  these  input  parameters  for  Do *  will  not  be  well-defined 
for  a  surface  comprised  of  random,  polydisperse  (but  monolithic)  aggregates,  the  average  distance 
between  aggregates  will  decrease  as  approximately  the  square  root  of  the  aggregate  concentration  (for 
monolayer-like  arrangements  of  aggregates  as  seen  in  Fig.  3),  while  the  aggregate  size  distribution 
remains  constant.  The  result  is  that  D0*  should  decrease  with  increased  aggregate  loading,  first  rapidly, 
and  then  gradually. 

By  examining  Fig.  1  more  closely,  it  can  be  seen  that,  once  the  Cassie  state  is  attained,  further 
increases  in  silica  loading  result  in  modest  increases  in  apparent  contact  angles.  In  a  composite 
interface,  the  higher  apparent  contact  angle  results  from  a  decrease  in  fa, total,  or  equivalently,  an  increase 
in  D *  (see  Supporting  Information  Section  SI  and  Refs  1-3  for  examples).  Thus,  the  experimental  data 
indicate  that  D*  increases  with  increased  FF-silica  loading,  contrary  to  the  expected  behavior  of  Do*. 
(A  more  quantitative  approach  using  equivalent  model  surfaces  is  presented  in  Supporting  Information 
Section  SI).  Therefore,  based  on  Eq.  (1),  u  must  increase  with  increased  FF-silica  content.  In  the 
specific  case  of  a  hierarchical  surface  with  well-defined  features  replicated  at  multiple  length  scales,  it 
has  been  proven  mathematically  that  u  will  increase  as  features  are  replicated  at  smaller  and  smaller 
length  scales.  The  fractal  dimensionality  of  the  textured  surface  also  represents  another  useful  single 
geometric  parameter  that  indicates  the  level  of  fine-scale  roughness  (see  Supporting  Information  Section 
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S4  for  an  explanation).  Therefore,  the  increasing  fractal  dimensionality  of  the  surfaces  with  increased 
FF-silica  content  seen  in  Fig.  5  implies  that  u  also  increases  with  increasing  FF-silica  content,  and  thus 
explains  the  observed  contact  angle  behavior.  An  increasing  value  of  u  should  also  stabilize  the  Cassie 
state  for  lower  surface  tension  liquids,  and  result  in  the  behavior  reflected  in  Fig.  2. 

In  addition  to  the  geometric  analysis,  a  careful  examination  of  the  cross-sectional  images  seen  in 
Fig.  4  revealed  a  qualitative  change  in  the  nature  of  the  topography  as  the  FF-silica  content  was 
increased.  (Note  that  these  changes  were  seen  over  the  course  of  examining  almost  100  separate 
images,  with  Fig.  4  providing  representative  examples).  At  loading  levels  below  40  wt%  silica,  there 
appears  to  be  an  “excess”  of  fluoropolymer  binder,  which  fdls  in  small  gaps  between  aggregates  and 
sub-aggregates  and  provides  a  relatively  thick  coating  of  the  substrate  that  partly  submerges  the  smaller 
aggregates.  The  result  is  a  surface  that  is  conformal  only  above  a  length  scale  of  a  micron  or  so.  Above 
40  wt%  FF-silica,  however,  there  appears  to  be  almost  no  “excess”  binder.  In  this  case,  the  highly 
confonnal  coating  preserves  the  roughness  of  the  aggregates  down  to  sub-micron  scales.  Although  the 
most  prominent  change  occurs  at  approx.  40  wt%  FF-silica,  where  the  aggregates  also  begin  to  clump 
together,  the  coating  can  be  seen  to  confonn  to  progressively  smaller-sized  features  as  the  FF-silica 
loading  increases.  The  gradual  spread  of  conformality  to  smaller  length  scales  due  to  an  increasingly 
severe  shortage  of  binder  thus  appears  to  be  the  mechanism  responsible  for  the  observed  changes  in 
fractal  dimensionality,  contact  angle,  and  minimum  surface  tension  of  contacting  fluids  exhibiting  the 
Cassie  state,  as  a  function  of  FF-silica  loading. 

The  above  discussion  illustrates  that  innovative  methods  can  be  utilized  to  systematically  alter 
the  topography  of  randomly  textured  surfaces  to  obtain  desired  liquid  repellence  characteristics.  Even 
in  the  context  of  spatially-variable  robustness,  which  appears  to  be  characteristic  of  such  surfaces  (see 
the  discussion  of  immersion  depth  and  evaporating  droplets  in  Supporting  Information  Section  S2),  it 
appears  that  the  contact  and  sliding  angles  can  be  altered  in  systematic  fashion  by  controlling 
parameters,  such  as  the  fractal  dimensionality  or  the  distribution  of  roughness,  across  multiple  length 
scales.  For  sprayed  on  surfaces  consisting  of  rigid  aggregates  with  a  polymeric  binder,  the  aggregate  to 
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binder  ratio  influences  these  parameters  due  to  a  mechanism  by  which  local  accumulations  of  binder 
mask  the  fine-scale  roughness  of  the  aggregates.  The  aforementioned  mechanism  is  expected  to  apply 
to  all  types  of  coatings  containing  rigid  aggregates  dispersed  in  a  soft  binder,  provided  that  significant 
migration  of  binder  into  large  interstices  is  avoided.  Because  the  mechanism  involves  a  simple  excess 
or  shortage  of  binder,  it  is  expected  to  be  insensitive  to  the  details  of  the  coating  process  itself,  as  long 
as  the  aggregates  are  deposited  in  a  monolayer-like  fashion. 

4.  Conclusions 

Coatings  comprising  a  sprayed-on  mixture  of  fluoroalkyl-functional  precipitated  silica  and  a 
fluoropolymer  binder  have  proven  to  be  an  especially  useful  class  of  liquid  repellent  coating.  The 
combination  of  low  surface  energy  materials,  aggregate  porosity,  and  rapid  deposition  via  spray  coating 
appears  to  prevent  the  build-up  of  binder  in  interstitial  regions,  preserving  re-entrant  curvature  across 
multiple  length  scales  and  thereby  enabling  a  wide  range  of  liquid  repellency,  including 
superoleophobicity,  to  be  attained  via  a  simple  and  scalable  process.  In  addition,  rather  than 
accumulating  in  the  interstices,  the  binder  becomes  widely  distributed  across  the  surface  of  the 
aggregates,  enabling  a  mechanism  in  which  a  simple  shortage  or  excess  of  binder  controls  the  extent  of 
coating  roughness  at  very  small  length  scales,  thereby  controlling  the  extent  of  liquid  repellence.  The 
exploitation  of  these  advantages  is  expected  to  lead  to  significant  progress  towards  large-scale 
commercialization  of  superoleophobic  surfaces. 

Supporting  Information  Available 

Properties  of  fhroroalkylsilane-functionalized  silica  (FF-silica,  Table  SI),  complete  contact  angle  data 
with  uncertainties  (Table  S2),  original  SEM  micrographs  of  surfaces  (Figures  SI  and  S2),  regular 
geometric  model  for  estimating  liquid  repellence  parameters  of  FF-silica  composites  (Section  SI), 
robustness  properties  of  FF-silica  composites  with  accompanying  video  (Section  S2),  derivation  of  the 
dependence  of  average  distance  between  particles  on  particle  loading  (Section  S3),  and  derivation  of  the 
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linkage  between  fractal  dimensionality  and  extent  of  fine-scale  roughness  (Section  S4).  This  material  is 
available  free  of  charge  via  the  Internet  at  http://pubs.acs.org. 
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Table  SI.  Physical  and  Geometric  properties  of  fluoroalkyl-functionalized  silica  aggregates  used 

o  1 

in  this  work. 


Silica  Type 

precipitated  (HiSil  233) 

Average  Diameter  (nm) 

22 

BET  Surface  Area  (m2/g) 

92 

BET  C  Constant 

21 

Water  Vapor  Uptake  (wt%) 

2.8 

Wt  %  Fluorine 

9.9 

Grafting  Density  (chains  nm'  ) 

1.6 

Graft  Layer  Molar  Volume  (cc) 

311 

Average  Thickness  of  Graft  (nm) 

0.8 

Approved  for  public  release,  distribution  unlimited. 


Table  S2.  Summary  of  apparent  advancing  and  receding  contact  angles  of  water  (y  =  72.1 
inN/m),  diiodomethane  (y  =  50.8  mN/m),  ethylene  glycol  (y  =  47.7  mN/m),  dimethyl  sulfoxide  (y 
=  44  mN/m),  rapeseed  oil  (y  =  35.5  mN/m),  hexadecane  (y  =  27.5  mN/m),  dodecane  (y  =  25.3 
mN/m),  and  decane  (y  =  23.8  mN/m)  on  silicon  wafers  sprayed  with  FF-silica  filled 
fluoropolymer  coatings. 


Water 

Diiodomethane 

Ethylene  Glycol 

Dimethyl  Sulfoxide 

Sample 

®adv 

®rec 

®adv 

^rec 

^adv 

^rec 

®adv 

®rec 

Viton 

125  ±  1 

83  ±2 

102  ±  1 

66  ±4 

102  ±  1 

56  ±6 

94  ±2 

53  ±4 

10  wt% 

141  ±2 

78  ±5 

1 14  ±  4 

52  ±5 

113  ±  5 

27  ±5 

108  ±3 

27  ±4 

20  wt% 

156  ±2 

151  ±3 

138  ±5 

35  ±6 

147  ±3 

25  ±6 

137  ±3 

13  ±5 

30  wt% 

160  ±2 

160  ±2 

146  ±4 

95  ±8 

153  ±5 

139  ±8 

145  ±5 

5  ±  5 

40  wt% 

159  ±  1 

158  ±4 

148  ±4 

133  ±6 

156  ±  2 

146  ±5 

156  ±2 

132  ±3 

50  wt% 

162  ±2 

162  ±3 

145  ±3 

136  ±8 

154  ±  4 

145  ±  13 

152  ±2 

133  ±9 

60  wt% 

161  ±2 

162  ±2 

152  ±2 

146  ±4 

155  ±  2 

145  ±7 

155  ±2 

139  ±  10 

70  wt% 

168  ±3 

169  ±2 

151  ±3 

148  ±4 

156  ±  4 

149  ±4 

157  ±  1 

148  ±4 

80  wt% 

165  ±1 

165  ±  1 

152  ±3 

151  ±1 

157  ±  2 

148  ±3 

160  ±3 

153  ±6 

Rapeseed 

Oil 

Hexadecane 

Dodecane 

Decane 

Sample 

®adv 

®rec 

®adv 

®rec 

^adv 

^rec 

^adv 

®rec 

Viton 

86  ±2 

42  ±3 

72  ±2 

39  ±2 

70  ±3 

37  ±4 

58  ±2 

24  ±3 

10  wt% 

95  ±3 

26  ±3 

82  ±3 

5  ±  5 

20  wt% 

120  ±3 

5  ±  5 

101  ±3 

5  ±  5 

30  wt% 

141  ±5 

5  ±  5 

126  ±3 

5  ±  5 

40  wt% 

154  ±5 

5  ±  5 

135  ±  4 

5  ±5 

50  wt% 

151  ±6 

136  ±  5 

144  ±4 

5  ±  5 

60  wt% 

154  ±2 

138  ±7 

155  ±  1 

131  ±4 

154  ±2 

89  ±  16 

70  wt% 

156  ±2 

142  ±3 

156  ±2 

131  ±5 

156  ±3 

147  ±6 

80  wt% 

156  ±2 

153  ±2 

155  ±3 

146  ±4 

162  ±2 

155  ±2 

77  ±3 

5  ±  5 

2 


Figure  SI.  SEM  micrographs  of  FF-silica  /  fluoropolymer  coatings  (“top  down”  view  with 
substrate  parallel  to  the  plane  of  the  micrograph).  The  images  in  Figure  3  in  the  main  text  are 
cropped  and  re-sized  so  as  to  make  for  identical  magnification  and  to  fit  the  available  space 
while  maximizing  viewable  details.  The  flattened  and  smoothed  features  at  20  wt%  FF-silica  are 
due  to  the  presence  of  a  large  excess  of  binder,  which  is  absent  at  40  wt%,  60  wt%,  and  80  wt% 
FF-silica  loading  levels. 
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Figure  S2.  SEM  micro-graphs  prepared  from  FF-silica  /  fluoropolymer  coatings  mounted  in 


epoxy  and  cross-sectioned.  The  images  in  Figure  3  in  the  main  text  are  rotated,  cropped,  and  re¬ 


sized  so  as  to  create  identical  magnification  and  orientation,  and  to  fit  the  available  space  while 


maximizing  viewable  details. 
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Section  SI.  Idealized  Model  to  Describe  the  Wetting  Behavior  of  Spray-Deposited  FF-silica 
Filled  Fluoropolymer  Coatings 


o  o  o 
o  o  o  o 

Q  O  O 

■< - - >■ 

I  |  I 

2  R  2D 


Figure  S3.  Left:  SEM  micrographs  of  spray-deposited  FF-silica  filled  fluoropolymer  coatings. 
Right:  Idealized  geometric  model  used  to  describe  the  surface  of  spray-deposited  FF-silica  filled 
fluoropolymer  coatings. 


Spray-deposited  coatings  of  fluoropolymer  filled  with  FF-silica  described  in  the  main  text 
demonstrate  increased  super-liquid  repellency  with  increasing  FF-silica  content  for  liquids  with 
decreasing  surface  tensions,  embodied  by  high  apparent  contact  angles  (>150°)  with  low  contact 
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angle  hysteresis  (<  10°)  and  corresponding  low  sliding  angle.  The  wetting  behavior  of  these 
surfaces  implies  the  formation  of  a  composite  solid-liquid-air  interface  between  probing  droplets 
and  test  substrates,  with  pockets  of  air  being  present  between  micron  and  sub-micron  surface 
features  that  possess  inherent  re-entrant  curvature.  This  system  can  be  described  by  the  Cassie- 
Baxter  relationship  in  the  following  form: 

cos<9  =  cos  Oe+^s -1  (SI) 

where  cost?/,  is  the  contact  angle  on  an  identical  surface  absent  of  roughness  (Young’s  contact 
angle),  r </,  is  the  roughness  of  the  wetted  area,  <f>s  is  the  area  fraction  of  the  liquid-air  interface 
occluded  by  the  solid  texture,  and  is  the  fraction  of  the  solid  substrate  in  contact  with  the 
liquid. S2 

For  a  model  textured  surface  with  a  known  surface  energy,  geometric  parameters  can  be 
included  in  the  Cassie-Baxter  equation  to  accurately  predict  the  apparent  contact  angle  of  a 
millimeter-sized  droplet  of  a  given  liquid,  assuming  that  0E  is  satisfied  locally  on  micron  and 
sub-micron  features.  "  Although  the  coatings  described  in  the  present  work  are  devoid  of  any 
regularly  arranged  geometric  features,  as  depicted  in  the  SEM  micrographs  shown  in  Figures  3 
and  4  of  the  main  text  and  in  Figures  SI  through  S3,  an  idealized  geometric  model  may  be 
adopted  to  elucidate  the  increased  super-liquid  repellency  as  a  function  of  FF-silica  content. 
Features  on  FF-silica  filled  fluoropolymer  surfaces  can  be  approximately  described  as  a 
hexagonally-packed  array  of  spheres  of  diameter  2 R  with  a  spacing  distance  of  2D,  as  shown  in 
Figure  S3.  Using  this  model  the  Cassie-Baxter  equation  takes  the  following  form: 

cos£*  =  -l+l/D*[7t/2V3(l+cos6k)2]  (S2) 

where  D*  =[(R+D)/R]2,  =2(l+cos6k)/sin26k  and  <f>s  =n  sin26,/./2V3D*.  From  this  relationship 
the  fraction  of  liquid  in  contact  with  surface  is  inversely  proportional  to  the  spacing  ratio,  D*. 
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The  effective  spacing  ratio  of  the  present  on  the  sprayed  surfaces  can  be  calculated  by  plotting 
cos  6 *  vs  [7t/(2V3)](1  -t-  cos  dE )2  and  fitting  a  straight  line  to  the  data  (shown  in  Figure  S4). 
The  slope  from  the  linear  fit  can  then  be  used  to  determine  the  spacing  ratio  between  spherical 
surface  features  reported  in  Table  S3. 
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Figure  S4.  Linear  fit  used  to  determine  D*  of  the  idealized  surface  having  the  characteristics  of 
the  FF-silica  /  fluoropolymer  coatings  (herein  denoted  as  “RCxx”  where  “xx”  denoted  the  weight 
fraction  of  FF-silica  in  the  coating. 
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Table  S3.  Calculated  spacing  ratios  (D  )  and  geometric  values  (D,  R)  for  spherical  features 
present  on  sprayed  coatings  of  fluoropolymer  filled  with  various  amounts  of  FF-silica. 


Sample 

D* 

D  (pm) 

/?(|im) 

30  wt% 

5  ±14 

360 

290 

40  wt% 

5.1  ±3.4 

320 

250 

50  wt% 

5.9  ±3.8 

250 

170 

60  wt% 

7.5  ±4.0 

210 

120 

70  wt% 

12.7  ±6.8 

100 

40 

80  wt% 

14.6  ±9.2 

90 

30 

The  values  in  Table  S3  are  essentially  the  design  parameters  for  a  surface  with 
hexagonally-packed  spherical  feature  with  identical  surface  chemistry  as  that  of  the  sprayed 
surfaces  reported  herein.  The  spacing  ratios  for  these  equivalent  surfaces  do  indeed  increase  with 
increasing  silica  content,  consistent  with  the  Cassie-Baxter  model.  However,  these  values  do  not 
correspond  to  any  obvious  features  or  spacing  ratios  present  in  SEM  micrographs  of  the  sprayed- 
surfaces,  or  to  any  trends  with  increasing  silica  loading.  In  fact,  from  top-down  and  cross- 
sectional  SEM  micrographs,  the  spacing  between  surface  features  decreases  with  increasing  FF- 
silica  content,  which  should  increase  the  wetted  fraction  when  adopting  the  model  presented. 
However,  as  mentioned  in  the  main  text,  the  presence  of  hierarchical  roughness,  with  multiple 
levels  of  re-entrant  curvature,  increases  the  amount  of  liquid  in  contact  with  micron-sized 
features  which  have  sub-micron,  essentially  non-wettable,  pores,  decreasing  the  overall  wetted 
fraction  with  increasing  FF-silica  loading.  The  hierarchical  nature  of  these  stochastic  surfaces 
prevents  a  direct  comparison  to  ideal  surfaces  modeled  to-date. 
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Section  S2:  Robustness  of  the  Partially  Wetted  State 

The  stability  of  a  given  composite  solid-liquid-air  interface  can  be  estimated  from  a 
dimensionless  “robustness  factor”  A  ,  which  is  essentially  a  ratio  of  the  pressure  required  to 
distort  the  liquid-air  interface  between  a  probing  liquid  and  test  substrate  to  induce  a  transition  to 
the  Wenzel-state,  and  the  minimum  pressure  across  the  interface  from  the  effects  of  gravity  and 
the  Laplace  pressure  within  the  droplet;  A  given  droplet  is  predicted  to  be  in  the  Wenzel  state  if 
A  <1  or  in  the  Cassie-state  if  A  >1.  From  this  theoretical  framework,  one  can  approximately 
predict  whether  a  surface  with  a  known,  regular  surface  geometry  will  be  able  to  support  a 
droplet  of  a  given  liquid  in  the  Cassie-state,  and  what  external  pressure  can  be  applied  before 
transitioning  the  droplet  to  the  fully-wetted,  Wenzel-state.  The  pressure  required  to  transition  a 
Cassie-state  droplet  to  the  Wenzel-state  can  be  measured  with  a  simple  experiment  taking 
advantage  of  the  Laplace  pressure  present  inside  of  a  droplet  of  liquid  (P Laplace  =  2y iJRdrop)-  As  a 
droplet  of  liquid  evaporates  on  a  surface,  the  diameter  of  the  droplet  will  correspondingly 
decrease,  increasing  the  pressure  imposed  on  the  liquid-air  interface.  The  pressure  imposed  on 
the  liquid-air  interface  between  the  substrate  and  droplet  will  increase  until  the  bulging  liquid 
meniscus  contacts  the  unwetted  solid  substrate  below,  transitioning  the  droplet  transitioning  the 
droplet  to  the  fully-wetted,  Wenzel-state.  Such  an  experiment  is  illustrated  in  Figure  S5. 
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Composite 


Fully-wetted 


Figure  S5.  Select  still  images  of  a  droplet  of  methanol  evaporating  on  a  superoleophobic  surface 
textured  with  micrometer-sized  features  possessing  re-entrant  curvature;  At  a  discrete  pressure, 
the  droplet  transitions  to  the  fully  wetted  Wenzel-state  (Copyright  2008  National  Academy  of 
Sciences,  USA).S3 

t=0  sec  I  t=5  sec  I  t=10sec|  lt=15sec| 


Figure  S6.  Select  images  of  a  chlorofonn  droplet  evaporating/infusing  into  an  80  wt%  FF-silica 
surface  displaying  the  progressive  collapse  of  the  metastable  solid-liquid-air  interface  as  a  result 
of  the  increasing  Laplace  pressure  within  the  droplet.  The  gradual  transition  from  the  Cassie- 
state  to  the  fully  wetted  Wenzel-state  is  accompanied  by  a  continuous  reduction  of  the  contact 
angle  and  an  apparent  transfer  of  air  from  within  surface  asperities  to  inside  the  droplet.  A  video 
has  been  provided  along  with  the  Supporting  Information  presented  herein. 
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Similar  experiments  were  conducted  in  an  attempt  to  quantify  the  super-liquid  repellency 
robustness  of  the  sprayed-coatings  reported  in  this  work.  However,  discrete  transitions  were  not 
observed  for  Cassie-state  droplets.  Instead,  a  gradual  decrease  in  contact  angle  was  seen  for 
Cassie  droplets  as  a  function  of  decreasing  droplet  size  (increasing  Laplace  pressure) 
corresponding  to  a  gradual  transition  to  the  wetted  state  as  smaller  and  smaller  levels  of 
hierarchical  features  are  wetted.  Representative  still  images  of  this  experiment  can  be  seen  in 
Figure  S7  for  a  droplet  of  chlorofonn  evaporating  on  a  spray-deposited  surface  of  fluoropolymer 
filled  with  80  wt%  FF-silica.  Another  interesting  observation  from  Figure  S6  is  the  formation  of 
air  bubbles  coming  out  of  surface  asperities  and  into  the  droplet  concomitantly  with  reducing 
droplet  size  and  contact  angle.  This  observation  provides  additional  evidence  that  the  transition 
to  the  Wenzel  state  does  not  occur  discretely,  as  with  model  textured  surfaces,  but  occurs 
gradually  for  the  sprayed  surfaces  suggesting  that  a  single  breakthrough  pressure,  or  A*  value, 
cannot  be  used  to  describe  the  surfaces  produced  in  this  work.  Attempts  to  measure  such 
breakthrough  pressures  more  directly  by  slowly  immersing  samples  in  increasing  depths  of  fluid 
also  produced  variable  results,  although  in  some  cases  stable  air  pockets  remained  after 
immersion  in  three  feet  of  water. 

Additionally,  the  formation  of  air  bubbles  evolving  from  the  occluded  surface  of  the 
substrate  suggests  that  the  pockets  of  air  under  the  droplet  are  effectively  trapped,  unable  to 
escape  parallel  to  the  surface  when  subjected  to  an  external  pressure.  This  observation  may  help 

to  explain  the  repellency  to  low  surface  tension  liquids  (such  as  hexadecane)  for  which  the 

•  •  •  9  m 
equilibrium  contact  angle  6  is  «90.  Although  the  surface  features  present  on  this  surface  do 

not  possess  ideal  geometry  for  non-wetting  surfaces,  the  liquid-air  interface  occluded  by  the 
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droplet  is  supported  by  the  increasing  pressure  of  the  isolated  air  pockets  as  the  sagging  of  the 
liquid- air  interface  increases.  Therefore,  the  future  design  of  textured  liquid-repellent  surfaces 
may  afford  increased  robustness  if  low-lying  spaces  between  surface  features  are  isolated  from 
one  another,  trapping  air  and  preventing  the  collapse  of  the  composite  interface  necessary  for 
super-liquid  repellency. 

Section  S3:  Simplified  Calculation  of  the  Dependence  of  the  Average  Distance  Between 
Particles  as  a  Function  of  Particle  Loading  in  a  Stochastic  Monolayer 

The  key  assumptions  underlying  the  calculation  are  as  follows: 

1.  At  the  largest  scale,  the  surface  texture  is  dominated  by  protrusions  caused  by  rigid 
aggregates.  Otherwise  the  surface  texture  sufficiently  far  away  from  the  aggregates  is 
relatively  flat  and  smooth  as  a  result  of  homogeneous  film  formation  of  the 
fluoropolymer  binder. 

2.  Interactions  between  nearby  particles  do  not  create  significant  new  surface  texture. 
Similarly,  the  surface  texture  does  not  depend  on  the  details  of  the  deposition  process. 
Thus,  the  number  of  protrusions  in  the  surface  is  directly  proportional  to  the  number 
density  of  particles. 

3.  The  overall  density  of  deposited  particles  is  low  enough  that  overlap  is  rare  and  can  be 
neglected.  In  other  words,  the  deposited  particles  form  a  sparse  monolayer. 

The  effect  of  particle  loading  is  computed  on  the  basis  of  solutions  with  a  constant 
concentration  of  solids  (co)  but  variable  proportion  (expressed  as  a  weight  fraction  wp)  of 
particles  to  total  solids  (particles  and  binder).  Though  polydisperse,  a  sufficient  number  of 
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particles  exist  such  that,  locally,  they  are  deposited  with  an  average  volume  (Vp)  that  changes 
insignificantly  from  region  to  region  within  the  deposited  coating.  The  micrographs  shown  in 
Fig.  3  of  the  main  text  indicate  that  the  above  conditions  are  satisfied  reasonably  well  for  the 
films  under  consideration.  At  higher  volume  fractions,  the  tendency  of  particles  to  form  small, 
overlapping  groups  does  begin  to  emerge,  but  for  the  present  analysis  this  effect  is  not 
considered. 

Based  on  the  variables  defined  above,  the  average  number  of  particles  deposited  per  unit 
area  of  the  surface  (Np)  is  given  by: 

Np  =  wf  Co  Vd/  Ad  pp  Vp  (S3) 

in  which  Vd  and  AD  represent  the  volume  of  solution  deposited  and  the  area  over  which  the 
deposition  occurs,  respectively,  and  in  which  pp  represents  the  true  density  of  the  aggregate, 
assuming  that  it  is  fully  wetted  by  solvent. 

If  a  plane  parallel  to  the  substrate  is  divided  on  the  basis  of  a  well-accepted  method  such 
as  Voronoi  tessellation  (that  is,  into  distinct  regions  defined  by  the  identity  of  the  co-ordinates  of 
the  nearest  particle  center  of  mass,  with  one  region  per  particle),  then,  on  the  basis  of 
assumptions  1  and  2,  there  should  be  exactly  one  large  protrusion  per  Voronoi  polygon  on  the 
surface.  The  distance  between  protrusions  can  then  be  defined  on  the  basis  of  the  characteristic 
length  ( ds )  of  these  polygons,  which,  given  that  they  are  two  dimensional,  is  simply  the  square 
root  of  their  mean  area  ( as ).  Since  the  entire  plane  is  covered  by  these  polygons,  the  number  of 
polygons  per  unit  area  is  the  same  of  the  number  of  particles  per  unit  area,  and  thus  the 
characteristic  area  of  each  polygon  is  simply  the  inverse  of  the  number  per  unit  area,  that  is 

as  =  \  / Np  thus  ds  =  (1  /  Np)  1/2  =  (AdPp  VplwfCo  VD)m  (S4) 
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For  rough  surfaces,  the  geometric  parameter  D*  used  in  predicting  liquid  repellency  is 
typically  defined  as  the  ratio  of  the  distance  between  feature  centers  (ds  in  this  case)  to  the 
feature  size  (in  one  dimension).  The  characteristic  size  of  each  feature  will  depend  on  the 
conformality  of  the  film  and  the  characteristic  geometry  of  the  aggregates,  but  in  general  should 
scale  with  the  quantity  Vp  .  Hence 

D*~(AdPp  Vp/wfCo  VD)m/ Vp1/3~  (Ad  pp  /  c0  Vd)1I2Vp1/6  w/1/2  ~  C,  w/m  (S5) 
where  C;  is  a  constant.  Thus,  under  the  assumptions  given,  D*  will  decrease  with  increasing  wp. 


Section  S4:  Simplified  Explanation  of  the  Relationship  between  Fractal  Dimensionality  and 
Fine-Scale  Roughness  on  a  Textured  Surface 


The  essence  of  the  box  counting  routine  developed  by  Shibuichi  'S4'  is  that  the  fractal 
dimensionality  (for  a  surface  formed  by  iterative  convolution  of  a  smaller  pattern  on  to  the  same 
pattern  repeated  at  a  larger  scale)  may  be  defined  as  1  +  r,  in  which  r  represents  the  ratio  of  the 
cross-sectional  path  length  found  by  examining  the  surface  on  the  smaller  scale  to  the  path  length 
traced  solely  by  the  feature  on  the  larger  scale.  A  simple  example  is  illustrated  below. 

If  a  fractal  surface  is  formed  by  taking  the  following  pattern 

ZA 


For  each 


substitute 


which  gives  rise  to  a  self-similar  surface  that  (up  to  3  iterations)  appears  as: 
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(Note  how  the  angle  of  the  “steepest”  surface  features  scales  with  the  number  of  iterations,  hence 
any  sufficiently  iterated  fractal  surface  should  have  re-entrance  with  a  minimum  re-entrant  angle 
approaching  0°).  In  the  above  example,  the  fractal  dimension  of  the  cross-section  (using  the 
Pythagoream  theorem  with  an  assumed  45°  angle  between  the  protruding  segments  and  the 
background)  is  4  /  (2  +  2  )  ~  4  /  3.4,  since  substituting  the  smaller  pattern  increases  the  path 
length  by  this  ratio.  The  Shibuichi  technique  (which,  as  stated  earlier,  estimates  the  fractal 
dimensionality  as  1  +  r)  would  then  yield  a  fractal  dimension  of  ~  1  +  4  /  3.4  ~  7.4  /  3.4  ~  2.2  for 
this  surface  (if  iterated  infinitely). 

The  Wenzel  roughness  of  the  surface  (area  ratio  analogous  to  r)  after  only  one  iteration 
would  also  be  ~  4  /  3.4,  assuming  the  surface  exhibits  unidirectional  fractal  topography  (that  is, 
the  surface  displays  a  uniform  cross-section;  for  surfaces  with  multi-directional  fractal 
topography  the  ratio  would  be  higher  but  would  follow  the  same  trend  outlined  below).  In  fact, 
the  Wenzel  roughness  is  given  simply  by  rn,  where  n  is  the  number  of  iterations,  in  the  case  of 
the  uniform  cross-section.  The  number  of  iterations  n  also  defines  a  convenient  scale  for  the 
roughness,  since  each  iteration  involves  shrinking  the  length  scale  (in  the  case  above,  by  a  factor 
of  ~  3.4).  It  is  clear  that,  for  any  value  of  r  modestly  greater  than  one,  the  function  rn  will  rapidly 
increase  as  n  increases,  or,  as  the  length  scale  decreases,  with  the  increases  becoming  more 
drastic  as  the  value  of  r  increases.  (A  similar  expression,  using  the  range  of  length  scales  rather 

cc 

than  number  of  iterations,  is  noted  by  Shibuichi  as  well  as  in  the  earlier  work  of  Onda.  )  Hence, 
for  any  surface  that  may  be  represented  reasonably  by  an  iterated  pattern,  it  is  reasonable  to 
conclude  that  as  the  fractal  dimensionality  (1  +  r)  increases,  the  roughness  present  at  small  scales 
will  also  increase. 
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Precipitated  silica  is  an  excellent  example  of  an  aggregate  that  can  be  represented  as  a 
fractal  system,  ’  and  the  surfaces  under  study  are,  to  a  varying  degree,  conformal  with  the 
silica.  Moreover,  these  surfaces  do  exhibit  a  characteristic  value  of  r  from  1-1.3  over  a  wide 
range  of  length  scales.  Hence,  so  long  as  n  is  reasonably  large,  an  increase  in  small  scale 
roughness  should  correlate  with  increased  fractal  dimensionality.  To  estimate  n  for  precipitated 
silica,  the  primary  particle  diameter  (22  nm,  from  Table  SI),  along  with  the  previously 
reporteds6,S7  fractal  dimensionality  of  1. 7-2.4  (estimated  here  as  2  for  simplicity)  were  used, 
along  with  the  characteristic  aggregate  size  (10  -  20  pm,  see  Figure  S2).  To  estimate  n,  the 
nature  of  the  iterated  structure  must  also  be  determined.  Since  both  the  primary  particles  and  the 
primary  aggregates  tend  to  be  spherical,  the  repeated  unit  should  be  a  spherical  sub-aggregate. 
Such  sub-aggregates  can  comprise  as  few  as  4-5  particles  or  sub-features,  but  likely  not  more 
than  10-20  based  on  the  lack  of  clearly  defined  sub-structural  levels  in  precipitated  silica 
aggregates.  If  No  is  the  number  of  sub-features  that  are  clustered  into  the  repeated  unit,  then  the 
total  number  of  particles  (Na)  in  the  complete  precipitated  silica  aggregate  is  given  by  Na  =  No". 
In  addition,  if  xo  represents  the  diameter  of  the  primary  particle  and  xa  the  aggregate  diameter, 
then  Na  =  ( xJxoY ,  where  d  is  the  fractal  dimensionality. 

Combining  and  taking  the  logarithm  yields 

n  log  N0  =  d{  log  xa  -  log  x0)  (S6) 

Using  base- 10  logarithms  and  the  values  mentioned  previously  yields  a  value  of  5-6  for  the  right 
hand  side  of  Eq.  S6.  For  No  =  4,  n  ~  8-10,  while  for  No  =  20,  n  ~  4-5.  Hence,  a  reasonable 
estimate  for  n  is  in  the  range  of  5-10,  which  is  certainly  large  enough  to  increase  the  roughness 
ratio  substantially  with  increasing  fractal  dimensionality.  For  instance,  if  r  =  1.2,  then  the 
Wenzel  roughness  rn  will  be  about  50%  greater  than  if  r  =  1.1  (2.5  vs.  1 .6)  when  n  =  5. 
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